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The properties of a class of topological quantum phase transition (TQPT) are analyzed based on 
a model proposed by Haldane. We study the effect of finite temperature on this phase transition. 
We have found that finite temperature would drive this TQPT to be a crossover, while it is stable 
against the weak short range interaction. When the interaction is strong enough, however, this 
TQPT is unstable and other states would emerge. Then we investigate the effect of the on-site 
energy in the original haldane model. The critical difference between our TQPT and the topological 
phase transition in conventional quantum Hall system is discussed. Finally, we discuss the potential 
application of our analysis to a topological phase transition proposed in a realistic system. 

PACS numbers: 71.10.Fd, 71.30.+h, 73.43.Nq 
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I. INTRODUCTION 

Landau's theories of the continuous phase transition 
have been successful and provided a paradigm in the 
condensed matter physics. The spontaneous symmetry 
breaking plays a central role in Landau's theories. In 
zero temperature however, it is possible that the change 
of fermi-surface topology or the quantum fluctuation in- 
duce new kind of phase transitions^ beyond Landau's 
paradigm. The different phases in these phase tran- 
sitions are not classified by different symmetries, in- 
stead, they are characterized by different topology of 
Fermi surface^, different quantum numbers^ 5 - or non- 
local topological order parameters' 7 ^ etc. This kind of 
phase transition, known as topological quantum phase 
transition(TQPT), has attracted considerable attention 
in recent years because they are closely related with the 
field of Topological Quantum Computation 9 . Despite of 
great efforts, many properties of the Topological phase 
and TQPT, such as its order and universality classes , 
the effect of temperature and interaction and its stabil- 
ity, are still far from being completely understood. For 
example, since there is no spontaneous symmetry break- 
ing in the TQPT, the Mcrmin- Wagner Theorem— - 1 could 
not be applied to preclude the possibility of a finite tem- 
perature phase transition in 2D system. In this paper, 
we have partly answered above questions. 



II. THE HALDANE MODEL 

In this paper, we present a theoretical analysis of 
a class of TQPT. Our starting point is the Haldane 
modelii, which was proposed to study the realization 
of "parity anomaly" in graphene-like condensed matter 
system and is the first work of anomalous Hall effect. 
The phase transition in this model, essentially, is the 
transition between two phases with different parity, or 
in another word, different Chern numbers, which is in- 
duced by adjusting parameter to change the sign of the 



effective mass of the " relativistic" electrons. We show 
that this TQPT is a third order quantum phase transi- 
tion, while most other conventional continuous quantum 
phase transitions are second order. We also analyze the 
effect of finite temperature on this TQPT, we find that 
there is no phase transition at T > 0, in another word, 
there is no finite temperature phase transition for the 
Anomalous Hall Effect (AHE). The temperature drives 
the TQPT to be a crossover. The critical behavior as 
well as the properties of correlation function in the finite 
temperature have also been analyzed, then an experi- 
ment is proposed to detect these unconventional proper- 
ties. When there exists short range interaction, we deal 
with this problem by the method of mean field and show 
that this TQPT is robust against the weak short range in- 
teraction for either the repulsive or attractive case. This 
is very different from the regular Fermi liquid. Strong 
interaction would change the pictures: strong enough re- 
pulsive interaction would lead to a charge density wave 
(CDW) statei 2 -, while strong attractive interaction would 
result in a superconductor 13 . To make our above results 
applicable to the orginal Haldane model, we investigate 
the effect of the on-site energy , which shifts the degen- 
eracy of the two fermi points. Furthermore, we point 
out the difference between this kind of TQPT and the 
disorder-induced localization-delocalization transition in 
the conventional quantum Hall systems^, as well as in 
the quantum spin Hall effec t 15 ' 16 . Finally, we show that 
though the Haldane model could not be realized in the 
present experiment, our results can be applied to a large 
class of TQPT which could be realized in experiment. 

First, we briefly review the Haldane mode l 11 ' 18 and 
focus on the phase transition in it. Haldane considered 
a (2+l)-dimensional ((2+l)-D) model in which spinless 
fermions hop on a half-filling honeycomb lattice and cou- 
ple with a periodic external magnetic field. There is no 
"net" magnetic flux through a unit cell and thus there 
is no Landau quantization. However, a next nearest- 
neighbor hopping term with a phase ^e 1 * is introduced in 
the honeycomb lattice, which is a time-reversal breaking 
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term and leads to a " chiral" fermions without the fermion 
doubling effect. This system provides a condensed- 
matter analog of (2+l)-dimensional electrodynamics 17 
due to its linear energy spectrum around the Fermi point. 

In the honeycomb lattice, noninteracting fermions be- 
have as a semimetal with the Fermi surface shrunk into 
to two incquivalcnt points in one unit cell as shown in 
Fig. 1. The general action of (2+l)-D noninteracting 
fermions in a uniform or periodic potential is given by 

r d 3 k 

S = Tr t^^WG-HW). (1) 

For the Haldane model, it is convenient to represent G _1 
in a basis of spinors ^>k — (^kAj^kB) 7 ^ with A and B 
representing the different sublattice of the honeycomb 
lattice. Therefore we have 
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FIG. 1: (a)The honeycomb lattice as a superposition of two 
sublattice A,B. The basis vectors are ai, aa, as, two sublattice 
are connected by bi,b2,b3. (b) The Brillouin zone, R,i , R.2 
are basis vector in reciprocal-lattice. K, K' are two distinct 
fermi point in one Brillouin zone 

G~ l = + 2t 2 cos $ ^2 cos ( k ' b i)] + d ' °"> ( 2 ) 

i 

where d • a = Y2i ^ cr * with d\ — t\ J^i cos(k • ai), d 2 = 
*i sin(k • ai) , d 3 = — i 2 sin $ ^ ■ sin(k • bi), and the 
Pauli matrices a l (i — 1,2,3) . Here we set the on-site 
energy M in the original Haldane model to be zero for 
convenience, we would discuss the effect of M in Sec- 
tion V. t\ and t 2 e 1 ^ are the tunnelling amplitudes be- 
tween the nearest and the next-nearest neighboring sites. 
The parameters ai , a2 and a3 are the displacement from 
one site to its nearest neighborhoods and bi = a2 a3, 
b2 = &3 — aij etc. We could calculate the Chern Num- 
ber (Hall conductance) for this model using the standard 
representation^ 

<Jxy = J d 2 k[n ■ d x ii x d y n] (3) 

with n = d/|d|. Finally, we could obtain 

(1 if t 2 sin $ > 
-1 if t 2 sin$ < . (4) 
otherwise 



The Chern number is the topological quantum number—. 
It was first introduced into the condensed matter physics 
in the integer quantum Hall effect to explain the stability 
of the Hall conductance to weak pertubations 20 . The 
physical meaning of the Chern number is the quantum 
Hall conductance in our case. Different Chern numbers 
characterize different topological phase and we focus on 
the phase transition between them. 

This phase transition could be interpreted in a tradi- 
tional way by calculating the singularly of the ground 
state energy.. The energy spectrum for (2) is : 

E(m, k) = ±^fdJ+dJ+4 (5) 

where di have been defined above, and m = t 2 sin(<I>). In 
the ground state, the lower band is filled completely so 
the ground state energy 

E{m) = J J d 2 k E(k,m) (6) 




is (c) - » - " - (d) 



FIG. 2: The ground state energy(a) and its first order deriva- 
tive^), second order derivative(c), third order derivative (d) 

To find the singularity of E(m) in the phase transition 
point m=0, we calculate E(m) and its derivative, the nu- 
merical result is shown in Fig. 2. We can see that the sec- 
ond deviative is continuous while the third one is not in 
the phase transition point, which means that this quan- 
tum phase transition is the third order quantum phase 
transition. This result could be explained in a heuristic 
way when we consider the behavior around the two dis- 
tinct independent Fermi points. These two independent 
effective Hamiltonians H a (a = 1,-1 represent different 
Dirac Fermi points) are given by 

H a = B a -a (7) 

where D Q = (akyC,k x c,m a ), c = 3/2ti, m a — 
— 3v3at2 sin$. The parameter c is the "light velocity" 
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in the relativistic Hamiltonian and below we set c 
for simplicity. 

The total Hamiltonian is^i : 



(8) 



^±1 represent fermions around different fermi points. 
When we adjust the parameter <E> from negative to pos- 
itive, the mass of both Dirac fermions change sign and 
the phase transition occurs. 

Because the two Dirac Fermi points are independent 
with each other, below we only focus on one of them. 
From the Hamiltonian (7), we can easily calculate the en- 
ergy spectrum -E'(k) = ±\/k 2 + m 2 with + for the "con- 
duction" band and — for the "valence" band. The energy 
of the ground state could be estimate as 



E(m) = - J J d 2 ^^/* 2 



-H 3 + C. (9) 



C is the physical cutoff of the integral and is obviously 
analytic in the phase transition point m = 0, because we 
only concern about the singularity of E(m) at the crit- 
ical point m = 0, C is unimportant. From (9) we find 
that the third order derivative of E(m) to m is discon- 
tinuous in the point m = 0, so this TQPT is a third 
order phase transition. The analysis above implies that 
Hamiltonian (8) provides a effective approximation of the 
actual Hamiltonian in eq. (1) if we only concern the be- 
havior around the Fermi point. Below we use eq.(8) to 
analyze the properties of the phase transition for simplic- 
ity. 



III. FINITE TEMPERATURE EFFECT 

A typical example of finite temperature effect on the 
quantum phase transition is the transverse Ising Model 1 , 
which is fermionized by the Jordon-Wigner transforma- 
tion. The finite temperature drive the quantum phase 
transition into a crossover. It is shown that there is a 
similar behavior in our case. When the temperature is 
not too high, the model could be expressed in terms of 
a continuum canonical fermion field W = ^ i) T ■ Its 
partition function is 



l/T 



Z F = J D* T (x,r)D*(x,r)ea;p(-y dr J d x.Lp) 

(10) 

with 



l f = ^i(^+a/ + 9/-^ ! )*-i (11) 
+* f i(^: + ~ V + ma*)*!. 

Here r is the imaginary time. From (10) and (11), we 
perform the standard scaling transformation by rescal- 
ing lengths, times, as well as the field etc . We define 







x' = xer l ,r' = Te~ zl ,~fy' — with e~ L being the di- 
mensional rescaling factor. To make sure the new Lp 
having the same form in terms of x' , t' , we find that 
z = 1 and m' — me 1 . So we have dim[m] = 1 and the 
critical exponent v = 1. 



(a)" 



(b) 



(d) 



FIG. 3: The free energy(a) and its first order derivative(b), 
second order derivative(c), third order derivative (d) at dif- 
ferent temperature 

At finite temperature T, the free energy F is given by 
d 2 k 



F = -T 



(27T) 



r 21n(cosh(\/k 2 + m 2 /T)). (12) 



To find out whether there is a finite temperature phase 
transition, we should calculate the derivative of the free 
energy. Numerical results are shown in Fig. 3. (Because 
we only consider the singularities around the phase tran- 
sition point, we set a physical cutoff of the integral and 
it would not effect the singularities of (12).) We can see 
that the free energy is analytic there is no singularity in 
(12) at T 7^ 0, thus no phase transition occurs. How- 
ever, there is a crossover just like in the transverse Ising 
model at finite temperature. The phase diagram is shown 
in Fig. 2. There are two regions in the phase diagram: 
the low-T massive Chiral-liquid and the continuum high- 
T region (quantum critical region). However, the tem- 
perature could not be too high that the continuum ac- 
tion Lp in (11) fails to provide a good approximation 
of the real system. Below we analyze the properties in 
both regions by calculating the fermion Green's function: 
G F (x, t) = (*(x, i)*(0, 0)). In the quantum field theory, 
the Green's Function for the massive Dirac Fermions at 
finite temperature is : 



Gf{ X , t) = {i^d^ + m) ab g F ( X , t) 



5F(x,t) 



d 2 k 



1 



-i\z.-x.-\-i\ / 'k. 2 -\-'rn' 2 t 



(2tt) 2 Vk 2 + m 2 1 + e^+^/T 



(13) 



(14) 
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We use the 2D representation of 7 matrix (/i=0,l,2) 
rather than the ordinary 4D representation, a, b is the 
spin of the Dirac electron, t denotes the real time. 

In our case, there are two kinds of Dirac fermions^i 
and $_i) and both of them contribute to the Green func- 
tion, the Green function is sightly different from the stan- 
dard formula (13), as we will show below. 



GiUtmin 



UidRmin 



FIG. 4: Finite -T phase diagram of the TQPT as a function 
of the mass m of the Dirac fermion and the temperature T 

Low-T massive Dirac-liquid, T <C \m\: In this region, 
the temperature is very low and thus the physics is con- 
trolled primarily by the critical line T = 0. First we 
analyze T = equal-time correlations given by 



0f(x,O) 



d 2 k 



(27T) 2 VF" 



dk kJo(k\x\) 
2^Vk 2 



m 
(15) 

We could approximately estimate the asymptotic behav- 
ior when |x| — > 00 by the contour integration which picks 
up contribution from the poles at which function is di- 
vergent in the complex k plane. For large x, those poles 
closest to the real axis provide the dominant contribu- 
tion to the integral. The leading result is G(|x|,0) « 
C(a;)e _ ' m " x ' (|x| — > 00), and thus the equal-time correla- 
tion exponentially decays in the spacial dimension with 
a correlation length £ = l/|m|. 

Then we calculate the equal-space correlation. For T = 
0, we have 



So G(0, t) decays with a power law and oscillatory 
with a frequency |m|. The correlation length in the real 
time axis (coherent time) is infinite when T = 0. We 
can see that though there is a gap, quantum systems at 
T = indeed have a long-range phase correlation in time 
which could not be seen when we map the D-dimcnsional 
quantum system to the D+l dimensional classic system 1 . 
When the temperature is very low, the density of the 
quasiparticles is small, so the life of a quasiparticle (co- 
herent time) holds for a long time in this region. 

Now we discuss the observable effect of G(0, t) in this 
region. We use the tunnelling effect in a bilayer Haldane 
model. There is a voltage V between the bilayer. By cal- 
culating the I-V curve in the our system, we can find that 
it is totally different from the ordinary electron systems 6 . 

Notice that we only consider the long-term behav- 
ior of G(0, t) in experiment, that is t ^ (actu- 
ally it is t ^> t~~t , we set h = 1 in our case). So 
\m\/t > l/t 2 , the term ^-e i|m|t is neglectable in (18), 
so G AA (0,t) cc M e iM* . The correlation of the tunnel- 
ing operator I(t) = C\ A (0, £)G2a(0, t) is : 



</(t)l+(0)) 



-G AA {t, 0)G\ 



AA (r 



t, 0) a 



t 2 



where 1, 2 denote the different layer and A, B are dif- 
ferent sublattices in each layer. Using the linear re- 
sponse formation, it is easy to find that the tunneling 
current / oc m 2 V, which is different from the massless 
case I cxV 3 . 



\HAf 



A2 



B2 



A2\ 



v 

T 



g(o,t) 



(2tt) 2 Vk 2 + m 2 



Doing this integration directly, we get 



1 



7(0, t) oc _ e <l"»l*. 
it 



So we find 



G AA {Q,t) cc 



it 2 



J\m\t 



J\m\t 



t 



(16) 



(17) 



(18) 



We notice that the term m in (13) doesn't appear in (18), 
because the m in different fermi points have opposite sign 
and cancel each other so G(0, t) is independent of the sign 
of m, which is reasonable in physics. 



FIG. 5: The bilayer Haldane model with voltage V between 
the layers. The I-V curve of the tunneling current could be 
measured in this device 

Continuum high T region, T> \m\: This region is also 
known as the " quantum-critical region" , which plays an 
important role in condensed matter physics because its 
potential relation with the heavy fermion systems 2 ^ and 
the high-Tc superconductivity^. The de Broglie wave- 
length of the excitation in this region is the same order 
as their mean spacing, so quantum and thermal effect en- 
tangle and both of them play an equally important role. 
We mainly focus on the line rn = 0, and calculate the 
equal-space correlation to see how the effect of tempera- 
ture change the behavior (13). From 
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7(0, f) = - 



d 2 k 



i\k\t 



(2tt) 2 \k\(l + e-\ k \/ T ) 



dk 



2^2 l + e -fc/T' 

This integral could be done by the contour integration. 
In the experimental condition: t 3> h/ksT (we set kg — 
H = 1 in our case) , it is easy to find that 



G(0,t) a T 2 e 



(20) 



We can see that (21) is totally different from (18). 
The equal-space correlation decays exponentially and the 
phase coherent time is proportional to 1/T in this region. 



IV. SHORT RANGE INTERACTION 

For fermions, the interaction and interaction-induced 
quantum phase transition play a central role in condensed 
matter physic o 26 ' 27 . We will show that our TQPT is sta- 
ble against both the weak repulsive and attractive inter- 
action, while strong enough repulsive/attractive interac- 
tion would lead to the CDW/superconductor state. 

First we consider the condition in the critical point of 
the TQPT (to — 0). The Hamiltonian for the spinlcss 
fermions in the honeycomb lattice with nearest-neighbor 
interaction V is given by 

H = £$>~(aj& i+a + h.c) + V{nl - \){n b l+SL \% 

i a 

(21) 

where a is the displacement from one site to its nearest 
neighbor site. For the repulsive interaction, we use the 
standard mean-field to deal with the interaction. Assum- 
ing {m) = 1/2 + l/2(-l) i A with i = for i G A and 
i = 1 for i E B, we get 

H = J2 £[-?( a &+» + h.c) + V{^- - A(-ln>, 



(22) 

Following the standard procedure^ 7 -, we calculate the en- 
ergy of the ground state. Minimize it with respect to A, 
we obtain the self-consistent equation 



1 = V 



d 2 k 



1 



(2tt) 2 ^/E 2 {k) + A 2 V 2 



(23) 



with E 2 (k) = Acos 2 {k x /2) + 4 : cos(V3ky/2ky)cos{k x /2) + 
1 at the critical point (to = 0). When we deviate 
from the critical point, there is a mass term in E 2 Qs), 
and it is estimated that E 2 (k, to) Acos 2 (k x /2) + 
4cos(v / 3fc y /2) cos(k x /2) + 1 + to 2 . 

By solving this equation, we find that there is a SM- 
CDW transition at a finite value V c (m) — 0.98. The 
numerical result is shown in Fig. 6. For the weak in- 
teraction V < V c , A = 0, and the TQPT is stable. 
It is consistent with previous study of the chiral sym- 
metry breaking induced by Coulomb interaction in Lay- 
ered Graphite 2 -. For attractive interaction, a mean field 



similar to the BCS theory could be done. The SM- 
Superconductor phase transition in the honeycomb lat- 
tice has been analyze d 13 ' 24 . These works support the 
result that, for the half-filling fermions, there exist a SM- 
Superconductor phase transition in V = — V c and when 
V > —V c , the SM is stable. This result is also correct 
when to deviates from slightly, so our TQPT is stable 
for the weak attractive interaction. 

In addition, it is interesting to study the effect of the 
temperature in the global phase diagram in Fig. 6: In the 
CDW phase, it is known that there is a finite-temperature 
phase transition for the 2D CDW system. When the tem- 
perature is high enough, the CDW order is destroyed. 
However, for the topological phase, as analyzed above, 
there is no thermodynamic phase transition at finite tem- 
perature, this result has shown a critical difference be- 
tween the topological phase and non-topological phase. 




FIG. 6: The phase diagram of the TQPT with the next- 
nearest interaction V, the solid line is the phase boundary 
between the CDW and topolqiogical phases, and the dot line 
is the phase boundary of our TQPT 



V. EFFECT OF THE ON-SITE ENERGY M 

Up to now, We assume the on-site energy in the orig- 
inal Haldane model M — 0, as a consequence, the two 
kind of Dirac Fermions ip\ and degenerate, the ef- 
fective mass of them: m a = — 3\/3at2 sin $. At the 
Phase transition point, the mass of both kinds of fermions 
change sign simultaneously thus the Chern number of the 
whole systme changes from —1 to +1. When the on-site 
energy M is included, as in the original Haldane model, 
the condition is different. Now the effective mass of the 
two kind of Dirac fermions m a = M — 3-\/3ai2 sin thus 
the energy of them no longer degenerate. The phase dia- 
gram with M has been demonstrated by Haldane (Fig. 7). 
We assume M > 0, as shown in Fig. 7, the original criti- 
cal point $ = is split into two different phase transition 
point $ Q = asin _1 (Af/3\/3i2), ct — ±1. At the point of 
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$+, for instance, only the mass of one kind of W_i change 
sign, while the other is not, thus the Chern number of the 
whole system is changed from +1 to 0, different from the 
condition when we set M — 0. At this phase transition 
point, we will find that only one kind of Dirac fermions 
experience the same singularity as shown in Fig. 2, while 
the energy of the other one is analytic in this transition 
point. However, if we only focus on the properties of the 
phase transition point, the analytic part is not impor- 
tant, so all the results in above sections can be applied 
to this kind of phase transition, that the it is a third 
order quantum phase transition, the temperature would 
drive it to be a crossover and the weak interaction can't 
change the TQPT. 




-4-3-2-10 1 2 3 4 



FIG. 7: The phase diagram in the original Haldane model. 
The on-site energy M shifts the degeneracy between the two 
fermi points and splits our original TQPT into two TQPTs 
(+-) 



VI. APPLICATION AND DISCUSSION 

Now we discuss the important difference between 
our TQPT and one of the most typical TQPT, i.e., 
the disorder-induced localization-delocalization transi- 
tion between different Hall plateaus in the Quantum Hall 
effec t 14 ' 29 . This kind of phase transition is controlled by 
the impurity and it is, in nature, a kind of percolation 



phase transition with a critical exponent v — 7/3. In our 
TQPT, however, there are no impurities and the phase 
transition is controlled by the mass of Dirac fermions. 
From the analysis above, we can get that v = 1. So our 
TQPT is totally different from the phase transition in 14 . 
It is also different from the localization-delocalization 
transitions in the quantum spin Hall effect 1 ^ with v = 
1.66. 

Since the Haldane model is an artificial model which 
could not be realized easily in experiments, we would turn 
to other physical systems in order to find this TQPT 
experimentally. In a recent work about quantum spin 
Hall effect£i2£, an experiment in HgTe Quantum well is 
proposed to realize this TQPT. This simplified model is 
based on the k • p perturbation theory. For each kind of 
electrons with spin | or J,, there is a Dirac- type subband 
due to the special structure of the quantum well. When 
the thickness of quantum well is adjusted to a certain 
point, the effective mass of two kinds of Dirac electrons 
would change its sign and a TQPT occurs. The system 
does not break the time-reversal invariance, thus there 
is no Hall conductance, instead, there exists a quantum 
jump of the spin Hall conductance Aa s xy = -jt- How- 
ever, there is a difference between this TQPT and our 
case. The mass term in£ is k dependent: m — m* — ck 2 , 
which makes their Hamiltonian different from the Stan- 
dard Dirac fcrmion. This difference results that the phase 
transition is no longer between two symmetric phase with 
the topological number a = ±1, but between a spin hall 
insulator (er s = 0) and a spin hall conductance (a s = 2). 
However, if we only focus on the phase transition, simple 
scaling analysis implies that the term k 2 is relevant in 
the phase transition point, so our results in this paper 
could be applied to this TQPT. 

It is natural to ask what happens if the interaction 
is long range (Coulomb interaction), which corresponds 
to the more physical situation. Since the weak short- 
range interaction is irrelevant in the critical point (m = 
0), simple scale analysis would imply that the TQPT is 
still stable for the weak long range interaction. However, 
when the interaction is strong enough, the situation is 
complex and whether there is a phase transition driven 
by the interaction is still controversial^!. To answer this 
question, it is necessary to do the renormalization group 
analysis of the (2+l)-Dimensional massive QED coupling 
with the gauge field. 
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